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Kibble-Zurek mechanism (KZM) uses critical scaling to predict density of topological defects and other 
excitations created in second order phase transitions. We point out that simply inserting asymptotic critical 
exponents deduced from the immediate vicinity of the critical point to obtain predictions can lead to results 
that are inconsistent with a more careful KZM analysis based on causality - on the comparison of the 
relaxation time of the order parameter with the "time distance" from the critical point. As a result, scaling of 
quench-generated excitations with quench rates can exhibit behavior that is locally (i.e., in the neighborhood 
of any given quench rate) well approximated by the power law, but with exponents that depend on that rate, 
and that are quite different from the naive prediction based on the critical exponents relevant for 
asymptotically long quench times. Kosterlitz-Thouless scaling (that governs e.g. Mott insulator to 
superfluid transition in the Bose-Hubbard model in one dimension) is investigated as an example of this 
phenomenon. 



The study of the dynamics of second-order phase transitions started with the observation by Kibble 1,2 that, in 
the cosmological setting, as a result of relativistic causality, distinct domains of the nascent Universe will 
choose different broken symmetry vacua. Their incompatibility, characterized by the relevant homotopy 
group, will typically lead to topological defects that may have observable consequences. 

In condensed matter (where the relativistic casual horizon is no longer a useful constraint) one can nevertheless 
define 3 5 a sonic horizon that plays a similar role. The usual approach to estimating the size of the sonic horizon 
relies on the scaling of the relaxation time and of the healing length that are summed up by the critical exponents. 
Critical exponents define the universality class of the transition, and this usually enables prediction of the scaling 
exponent that governs the number of the generated excitations (e.g., the density of topological defects) as a 
function of the quench timescale Tq for a wide range of quench rates. 

Here we point out that this simple procedure fails in an interesting and unexpected manner for the Kosterlitz- 
Thouless universality class. That is, one can expect that - in the asymptotic regime where the transition is 
extremely slow - critical exponents will suffice for such predictions. However, while for the quench rates attain- 
able in the laboratory one may still expect an approximate power law that relates density of excitations to the 
quench rate, the exponent that characterizes it will begin to approach predictions based on the critical exponents 
only asymptotically, and for unrealistically (one might even say, astronomically) large values of the "sonic 
horizon". Nevertheless, we show that the application of KZM can lead to predictions that are valid before the 
asymptotic regime characterized by the critical exponents becomes relevant. 

Timescale t at which the "reflexes" of the order parameter of the system, quantified by the relaxation time z, are 
too slow for its state to remain in approximate equilibrium with its momentary Hamiltonian (or free energy) 
controlled from the outside by the experimenter plays a key role. It is obtained from the equation 3-5 : 



:(e(i))=e(i)/e(i) 



(1) 



that compares relaxation time z with the rate of change of the dimensionless distance from the critical point, e.g. 
e=(T c — T) /T c where T c is the critical temperature. When e(r) is taken to vary on a quench timescale Zq as 
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equation (1) leads to: 



e(t) = t/z Q 



<*))=?. 



(2) 



(3) 



In phase transitions where the critical slowing down and critical 
opalescence can be characterized by power law dependencies of 
relaxation time and healing length, 



x(e) = x 0 /\er, C(e) = WM\ 
equation (1) can be imidiately solved: 

f = To(TQ/T 0 )ra% e=(z 0 /z Q )^ 



(4) 



(5) 



Above, v and z are the spatial and dynamical critical exponents that 
characterize the universality class of the transition, while t 0 and £ 0 
are dimensionful parameters determined by the microphysics. This 
leads to the characteristic scale 



£ = £o(t q /t 0 ) 1+ 



(6) 



It gives the size of the domains that break symmetry in unison, and, 
hence, dictates the density of topological defects left behind by the 
transition. 

Basic tenets of the above Kibble-Zurek mechanism have been 
confirmed by numerical simulations 618 , and, to a lesser degree 
(and with more caveats) by experiments 19 36 in a variety of settings. 
Refinements include phase transition in inhomogeneous systems 
(see 37 for recent overview) and applications of KZM that go beyond 
topological defect creation (see e.g. 3841 ). Recent reviews related to 
KZM are also available 42 - 46 . 

Our aim here is to note that when the critical slowing down is 
given by a more complicated dependence then the simple power law 
of Eq. (4), the resulting t and, therefore, | will vary in a way that 
cannot be fully characterized by the critical exponents that otherwise 
suffice to predict their scaling with the quench rate. That is, topo- 
logical defects or other excitations left behind by the quench will 
approach the scaling predicted by the critical exponents, Eq. (16), 
only asymptotically, and begin to conform with it only in the regime 
of extremely slow transitions that may be well out of the reach of 
laboratory experiments. In the regime of faster quenches that may be 
accessible to experiments a power law may still be locally a reasonable 
fit, although its exponent will vary slowly, approaching the asymp- 
totic prediction only very gradually. 

Results 

Kibble-Zurek mechanism in the Kosterlitz-Thouless universality 
class. This conclusion about the local power law dependence that 
approaches scaling dictated by the asymptotic vales of critical 
exponents is exemplified by the Kosterlitz-Thouless (KT) 
transition 47 49 . There the non-polynomial scaling of the healing 
length 

(7) 



e = e„exp [a 



where a ~ 1, is captured by stating that the spatial critical exponent 
v = °°, see e.g. 50 . This is a brief and dramatic way to sum up the faster 
than polynomial divergence of f , but it may tempt one to misuse Eqs. 
(15,16). Thus, formally, one could insert v = 00 relevant for the KT 
universality class into Eq. (16) to obtain: 

a = f 0 (T Q /T 0 ) 1/Z . (8) 

This equation may be (as we shall see below) asymptotically valid, but 
is unlikely to have the same range of validity asEqs.(15,16) regarded 
as the consequence of Eq. (1). In particular, for large Tq the exponent 

approaches 1/z, reflected in Eq. (8), only gradually. 



To see why, consider the equation for the relaxation time z oc £ 2 in 
the KT universality class: 

z(e) = t 0 exp (za j v /je| J (9) 
and assume, as before, 6= t/zQ. Equations (1) and (9) yield 

t 0 exp (^za j \Jt/z^j =i. (10) 

Thus, t is now a solution of a transcendental equation. It can be 
obtained as 



?/t 0 = (tq/to) 



(za/2) 



W[(za/2)^/z Q /z 0 _ 



(11) 



where W is the Lambert function. The above solution is plotted for 
different values of za/2 in Fig. 1A,B. This relation has been derived 
before and tested by numerical simulations in a 2D classical model in 
Ref. 51. 

Figure 1C shows that the slope of unity for the dependence of t on 
Tq (and, therefore, \ on Zq z ) is attained only for Tq many orders of 
magnitude larger than t 0 - for exceedingly slow quenches that are 
unlikely to be experimentally accessible. For even reasonably slow 
quenches the effective power law would be significantly less than 1, 
typically as small as —0.5 for tq ~ 10t 0 , gradually increasing to 
0.8. ..0.9 as Tq/t 0 grows to ~10 10 or so. 

Therefore, in transitions that exhibit KT-like non-polynomial 
scalings and result in symmetry breaking, the asymptotic behavior 
one might have inferred from the critical exponents sets in only in 
the regime that appears to be out of reach of experiments. For 
instance, the system would have to be large compared to the 

(f~10 10 ^o> which means (when we take modest Co ~ 10~ 10 m) 
that the size of the homogeneous system undergoing the transition 
should be large compared to say ~10 3 <f, or, in other words, 
kilometers! 

A similar difference between the critical limit and the critical 
regime, although with less dramatic consequences, arises near the 
para-to-ferro transition in the random Ising chain 53-55 : 

(12) 



H=-£K<jf +1 -5>«f, 



where /; and hi are randomly chosen ferromagnetic couplings and 
transverse magnetic fields respectively. Here in turn v = 2 is a solid 
number and it is the dynamical exponent that diverges in the critical 



regime 



2|d 



(13) 



The limit e— >0, where z — » °°, implies t, ~ Tq, i.e., a correlation length 
that does not depend on the quench time at all. However, a more 
careful analysis of the equation (1), employing the full formula (13) 
instead of just its critical limit, leads to a prediction that there is 
actually a slow logarithmic dependence on Tq, a conclusion that 
was confirmed by simulations in Refs. 53-55. 

A similar care proves beneficial for a non-linear quench 



e r = 



sign(r) 



(14) 



considered e.g. in Ref. 56. Here sign is the sign function and r > 0 is 
an exponent. Equation (1) yields 



l+zv 



"t t 0 (tq/t 0 ) 1+ ", e~ (t 0 /tq) 1+ », 
and the characteristic scale of length 



(15) 
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Figure 1 | In the textbook version of the Kibble-Zurek mechanism, the time t when the time evolution ceases to be adiabatic satisfies a power law 

' ,z '. In a log-log plot this power law becomes a linear function log 10 (f/To) = -j— — log 10 (TQ/To) + const, where t 0 is a characteristic timescale of 



the system. In (A), we plot i for a Kosterlitz-Thouless transition in function of Tq over many decades of the argument. This function may appear linear 
locally, i.e., in a range of one or two decades, but it actually becomes linear only for very slow quenches, and, consequently, for "astronomical" values of 
the frozen-out domain size J, Eq. (16). Indeed, in (B), we focus on the narrow range of Tq = 10° ■ 2 t 0 that are small enough for a realistic experiment. These 
plots maybe reasonably approximated by linear functions. In (C), a local slope dlog 10 (t/ To) /^l°gio( T Q/ T o) of the log-log plot in panel A in function of 
Tq. The slope 1, predicted in the critical limit when formally v — > ca , is achieved but only for Tq in the "astronomical" regime. When we focus on more 
realistic Tq, as in panel D, the local slope turns out to be significantly lower than in the critical limit. 



I^oOqAo) 1 



(16) 



Again, this simple but careful argument leads to the same conclusion 
as the calculations in Ref. 56. 

Kibble-Zurek mechanism in the Bose-Hubbard model. We 

emphasize that KT scaling is encountered in systems other than 
the classic KT transition in two dimensions (in which generation 
of vortex pairs occurs via thermal activation as the system is 
heated). Thus, while for the sake of definiteness, the discussion 
above was in the framework of finite temperature phase 
transitions, the universal character of the arguments makes the 
conclusions applicable also to quantum phase transitions in the 
ground state at zero temperature. The most celebrated example of 
the quantum KT universality class with z = 1 is the ID Bose- 
Hubbard model 57 : 



H= (bj +1 b, + bjb l+ ^ + ~u^n,[m-i) 



(17) 



where b\ is a bosonic annihilation operator at site I and wj = b\b\ is an 
occupation number operator. At a commensurate filling of 1 particle 
per site, the ground state of the model undergoes a K-T quantum 
phase transition from a localized Mott phase at / < J c to a superfluid 
phase at / > J c . The energy gap on the Mott side of the transition is 



A = A 0 exp 



(18) 



Here x = J/U is a dimensionless ratio of the hopping rate / to the 
interaction strength U. Using Ref. 58, it is possible to estimate: x c = 
0.26, A 0 = 0.2 /, and a = 0.3. 



Any quench from the Mott to the superfluid phase can be linear- 
ized near the phase transition 



e(t) = 



= -t/tQ. 



(19) 



The evolution ceases to be adiabatic at t = — t when the reaction time 
A -1 of the system equals the time remaining to the transition \ t\: 



Its solution is 



exp 



An / A = - 



= A n f. 



(tq/t 0 ) 



W[\Aq/t 0 J 
where the characteristic timescale is 

Ax, 



t 0 = 



a 2 A 0 ' 



(20) 



(21) 



(22) 



This inverse gap is proportional to the correlation length set at — i: 

Z~c;Ja 0 /a) Z = £ 0 (a 0 /a). (23) 



This correlation length is plotted in Figure 2. 

To summarize, the equation (7) applies in the critical regime 
where e<Cl and not only in the limit e— >0. When the last limit is 
taken in, say, the Bose-Hubbard model, then the equation implies a 
steep power law <f ~ T y , but a careful application of Eq. (7) in the 
whole critical regime shows that the steep power law is reached only 
for rather "astronomical" values of Tq and, especially, of C that can 
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Figure 2 | In (A), a log-log plot of the correlation length £, in function of the quench time Tq. In the textbook Kibble-Zurek mechanism there is a power 

, k V/(1+ZV 

law CQCTq 



In a log-log scale this power law would look like a linear function: log 10 [\ j <J 0 J = j^p — l°gio( T Q/ T o) + cons t- Our non-linear log-log 
plot can be reasonably approximated by a linear function locally, i.e., over a range of one or two orders of magnitude, but a local slope of this linearized 
approximation depends on the order of magnitude of Tq, as shown in panel (B). Fig. B shows the local slope dlog 10 ^C / £ 0 j / rflog 10 (TQ/To) ofthe log-log 
plot in panel A in function of Tq. For Tq — » °° the slope tends to 1, as predicted in the critical limit, but for any Tq that is reasonable experimentally it is 
significantly less than 1. For instance, the slope 0.9 is eventually reached at the "astronomical" Tq~ 10 10 To, but for a reasonable Tq = 10° ■• 2 t 0 the 
slope drops to a mere 0.2 ... 0.5. 



hardly be achieved in a realistic experiment. For more realistic 
quench times there is no power law, although in a narrow range of 
Tq there may appear to be one but with a much reduced exponent. 

Discussion 

We have seen that, in some cases, using KZM requires more than just 
inserting critical exponents (that are valid only asymptotically close 
to the critical point). Rather, to estimate the scale £ one must make 
sure that the key idea behind KZM - the scaling of the sonic horizon 
that results from the critical slowing down - is accurately described by 
the critical exponents in the regime probed by the experiment. This 
may seem like a straightforward requirement, but, as we have seen, 
there are situations where it may not be easy to satisfy. 

The example with Kosterlitz-Thouless scaling we have just dis- 
cussed may be extreme in that the scaling represented by the 
asymptotic values of critical exponents is attained only in the limit 
that is - FAPP - unreachable in the laboratory. Nevertheless, the 
KZM-like analysis based on the actual dependence of the gap on e 
enables prediction of the scaling modified to suit the range of the 
experimentally implementable quench rates. 

Key quantity for such considerations is e, the point where the 
behavior of the system changes character, and the corresponding t 
that defines the sonic horizon. However, even before one evaluates 
such subtleties exemplified by the KT transition, it is useful to verify 
the KZM prerequisite, i.e., whether transition starts and ends suffi- 
ciently far from the critical point to justify appeal to KZM. In experi- 
ments that involve emulation of condensed matter systems using e.g. 
trapped ions or BEC's and optical lattices this may be far from 
straightforward, as experimental constraints may force relatively 
short quench timescales (i.e., modest values of Tq/t 0 ) which means 



that e may be too large - sonic horizon may be defined too far from 
the critical point - to expect near-critical scalings to be relevant. 
Similar remark applies to sizes of systems: Unless sonic horizon 
~ £ is small compared to the size of the system, scalings predicted 
by homogeneous KZM will not apply (although - given certain addi- 
tional assumptions - one may be able to deduce their modified 
versions 37 ). 

A related and interesting issue is how does KZM fail when the 
assumptions are only approximately satisfied or even violated. 
Experiments such as 5 * suggest that this might be a "soft failure", 
i.e., some features of KZM (e.g., power law dependences) may still 
apply even while detailed predictions (exponents of these power 
laws) are unlikely to hold. 

There are also indications that even when the requirement of 
starting and ending the quench on the outside of the [ — e, + e] inter- 
val is satisfied only on one side, KZM like scaling may still emerge. 
While this is beyond the scope of the original KZM assumptions, it is 
clearly worthy of a more detailed investigation. 

Indeed, the Bose-Hubbard model is "gapless" on the superfluid 
side, so in this sense only the -eon the Mott insulator side is well 
defined. Yet, recent experiment suggests that 59 that power laws may 
approximate the post-quench state of the system, although (at vari- 
ance with KZM) their slopes appear to depend on where the system 
starts and ends the quench. Given that the investigated quench times 
were short (tq~to), so that quenches likely started and/or ended 
inside the [ — e, +e] interval, this is no surprise. 

One further complication that is worth noting is that the "ori- 
ginal" KZM was focused on predicting densities of topologically 
protected objects. More recent extensions use it to predict other 
properties of the system following continuous phase transitions. 
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Thus, the size of the sonic horizon has been used to estimate coher- 
ence length in the post-transition Bose-Hubbard superfluid. This is, 
again, an interesting extension, and there are settings (e.g., quantum 
Ising) where the numerical results (e.g., behavior of entanglement 
entropy 38 ) has been observed. However, as one moves away from the 
stable and well defined topological defects in integrable systems to 
less well defined and less stable characteristics (like coherence length 
in Bose-Hubbard systems that exhibit more complicated behavior), 
KZM may still yield useful "guidelines", but regarding it as prediction 
without further justification requires courage. 
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